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We investigate the growth of a film of some element Bon a substrate made of another substance A 
in a model of molecular beam epitaxy. A vertical exchange mechanism allows the ^4-atoms to stay 
on the growing surface with a certain probability. Using kinetic Monte Carlo simulations as well as 
scaling arguments, the incorporation of the A's into the growing B-layer is investigated. Moreover 
we develop a rate equation theory for this process. In the limit of perfect layer-by-layer growth, the 
density of yl-atoms decays in the B-film like (distance from the interface) -2 . The power law is cut 
off exponentially at a characteristic thickness of the interdiffusion zone that depends on the rate of 
exchange of a B-adatom with an A-atom in the surface and on the system size. Kinetic roughening 
changes the exponents. Then the thickness of the interdiffusion zone is determined by the diffusion 
length. 



PACS numbers: 68.35.Fx, 81.15.Aa, 81.15.Kk 
I. INTRODUCTION 

Heterolayers, where e.g. ferromagnets are in contact 
with antiferromagnets, semiconductors or superconduc- 
tors, give rise to new ordering and transport phenom- 
ena, which depend crucially on the interfacial struc- 
ture. Examples for such ordering phenomena are the 
exchange bias [1], or the cryptomagnetism [2]. Elec- 
tronic transport through a ferromagnetic-nonmagnctic- 
ferromagnetic sandwich ("spin valve" geometry) gives 
rise to the giant magneto-resistance [3,4]. Another 
example is the recently predicted possibility to en- 
hance or reduce a Josephson current magnetically by re- 
placing the tunnel barrier of a Josephson junction by 
a ferromagnetic-insulating- ferromagnetic sandwich [5]. 
These phenomena belong to the growing field of spin- 
tronics [6], where the spin degree of freedom is used for 
electronic signal processing. Interfacial mixing affects all 
of them [7]. For example, it leads to spin scattering dis- 
turbing the spin dependent transport properties. 

Therefore it is important to be able to control the var- 
ious physical processes that may spoil well defined inter- 
faces. Some of them proceed after growth such as bulk 
interdiffusion or chemical interface reactions like silicidc 
formation. However, there are also important processes 
taking place exclusively at the surface: For example the 
substrate may partially behave like a surfactant, when 
one grows a different material on it. It is this latter mech- 
anism which we investigate in this paper. The questions 
we want to answer concern the asymptotic concentration 
profile, the width of the interdiffusion zone and possible 
correlations among the substrate impurities within the 
growing layer. 

Specifically we consider growing some material B on 
a substrate A. Obviously the interfacial mixing requires 
that some substrate (A-) atoms get replaced by 5-atoms 
and "float up" on the surface until they get incorporated 
into the growing film. There is ample experimental evi- 
dence that such a behaviour occurs in very different sys- 



tems like Cr on Fe [8], AlAs on GaAs [9], Nb on Fe [10] 
or Au on Fe [11]. This process depends on several impor- 
tant parameters including the lattice mismatch and the 
interaction between the different atoms including mag- 
netic contributions. In particular the explanation of any 
ordering of A- and 5-atoms close to the interface would 
require a detailed investigation of these interactions [12]. 

The situation becomes considerably simpler, however, 
if one is interested in the physical properties further away 
from the interface. Then the concentration of ^4-atoms 
may be regarded as sufficiently low that their interaction 
as well as ^4-iJ-ordering become unimportant. The focus 
on this region justifies our simplified model, in which 
the interdiffusion zone depends only on the deposition 
rate F of _B-atoms, the diffusion constants Da and Db 
of the adatoms of type j4or B, respectively, and the rate 
E for the exchange of 5-adatoms with 4-atoms. The 
limit E I Db — > oo, where a 5-adatom exchanges with 
the first ^4-atom it encounters, would be realized, if the 
B-adatoms diffuse by an exchange mechanism [13], while 
the yl-adatoms diffuse by hopping. In the present paper 
we simplify the model even further by assuming that both 
kinds of atoms diffuse equally fast on the surface, Da = 
Db = D, with a diffusion constant independent of the 
surface composition. 

Apart from the exchange there is a second crucial in- 
gredient in the model: The ^4-atoms behave only partially 
as a surfactant in the sense that they can be overgrown 
by island edges. By contrast a perfect surfactant atom 
should "float up" also in front of an advancing island 
edge. 

Naively one would expect an exponential decay of the 
density profile of ^4-atoms far from the interface. It is 
the main result of this investigation that this is not al- 
ways the case: The incorporation of yl-atoms is much 
slower, giving rise to a power law decay of the concentra- 
tion profile in the limit of perfect layer-by-layer growth, 
D/F — > oo. In this case the width of the interdiffusion 
zone diverges, provided there are no finite size effects. By 
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contrast, we shall show that for finite D/F the width of 
the interdiffusion zone is no longer infinite, but a power 
law of D/F. 

This paper is organized as follows. In the next sec- 
tion we are going to define a simple solid-on-solid (SOS) 
model for epitaxial growth of a 5-layer on an yl-substrate, 
which allows for the irreversible exchange of 5-atoms 
with ^4's at the surface. In this model the A- atoms on 
the surface turn out to cluster in a time-periodic self- 
organized way, which is explained in Sec. III. The next 
three sections, Sec. IV - VI, are devoted to the limit 
D/F — > oo. First, in Sec. IV, we present a simple mean 
field argument leading to the prediction, that the concen- 
tration of ^4's decays algebraically in the 5-laycr. For a fi- 
nite system size this power- law is cut off leading to a finite 
width H of the interdiffusion zone, which is discussed in 
Sec.V, where also a scaling ansatz for the surface concen- 
tration ca is proposed. This scaling ansatz is confirmed 
by simulation results for one- and two-dimensional sur- 
faces in Sec. VI. The remaining sections deal with inter- 
diffusion for finite D/F. Sec VII contains simulation re- 
sults and scaling arguments, and in Sec. VIII a rate equa- 
tion theory is developed for the interdiffusion problem. 
In the Appendix we describe a very efficient implementa- 
tion of the simulation model for one-dimensional surfaces 
in the limit D/F — > oo. 

II. THE MODEL 

In order to model heteroepitaxial growth of a 5-layer 
on an ^4-substrate including the interfacial mixing, we 
introduce a simple solid-on-solid (SOS) model, where 
the lattice mismatch and most interactions between the 
atoms are neglected. Moreover we describe the exchange 
mechanism of 5-adatoms with ^4-atoms in the surface 
simply by a phenomenological constant exchange rate E, 
although it depends in reality on the system parameters 
as well as the local environment. 

The model is defined on a simple cubic lattice by the 
following kinetic rules (cf. Fig. 1): 

(1) Starting from an initially flat substrate consisting of 
^4-atoms, 5-atoms are deposited at randomly selected 
sites on the surface with deposition rate F. 

(2) As long as they do not have a lateral neighbor, the 
5-atoms diffuse on the surface with diffusion constant D. 

(3) When such a 5-atoms happens to sit on top of a A- 
atom, it can exchange vertically with rate E or continue 
to diffuse with rate D. 

(4) After an exchange, the B-atom stays irreversibly 
bound, whereas the ^4-atom diffuses on the surface with 
diffusion constant D. 

(5) There is no back exchange, when a ^4-atom sits on 
top of a 5-atom. 

(6) When two adatoms, regardless of their type, meet, 
they form a stable non-moving nucleation center of an 
island. 



(7) When an adatom, regardless of its type, reaches a site 
adjacent to an island, it is irreversibly bound, increasing 
the size of the island. 

(8) Both types of particles can diffuse down across terrace 
edges without being hindered by an Ehrlich-Schwoebel 
barrier. 

(9) There are no overhangs, i.e., we assume SOS growth. 

(10) The exchange of a 5-atom with an ^4-atom under- 
neath is forbidden, if the 5-atom is already part of an is- 
land, i.e., if it has a nearest neighbor at the same height. 
Hence yl-atoms can be overgrown by island edges. 

Measuring time in units of monolayers (ML) and 
lengths in units of the lattice constant a, this model is 
controlled by the two dimensionless parameters D/Fa 4 
and te = a 2 E/D. (In the following we set a = 1.) 
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FIG. 1. The growth process of the model. Batoms are de- 
posited with a deposition rate F onto the substrate. Adatoms dif- 
fuse on the surface and down to the binding site at an island edge 
with diffusion constant D irrespective of their type. In the top 
layer mobile B's can exchange with ^4's with an effective exchange 
rate E. 

This model will be investigated for one- and two- 
dimensional (d = 1 and d = 2) surfaces in the fol- 
lowing. Note that it reduces to the usual model for 
MBE growth (for a recent review sec e.g. [14]), if one 
does not distinguish the two particle types. Therefore, 
the D/F-dependence of quantities like island density, 
adatom density, surface width, etc. are the same as 
usual. For general values of D/F and E/Dwe use kinetic 
Monte Carlo simulations [15,16] in order to investigate 
the model. However, for one-dimensional surfaces in the 
limit D/F — > oo we implemented a much more efficient 
algorithm, which is described in the Appendix. 

III. CORRELATIONS OF THE A-ATOMS 

One of the most intriguing qualitative properties of this 
model is the time-periodic self organization of ^4-clusters 
on the growing surface with a period of one monolayer. 
Fig. 2 shows that the ^4-atoms (black) are first clustered 
around the nucleation sites of a new layer, but migrate 
towards the holes remaining in that layer, when the is- 
lands coalesce. Thus the characteristic distance between 
these clusters agrees with the typical distance between 
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the nucleation sites, the diffusion length 

£ D ~ {D/F)\ (1) 

as long as layer-by-layer growth persits. This can be 
verified by examining the lateral correlations of A-atoms 
on the surface after deposition of t monolayers: 

1 Ld 

g(f,t) = Yd^2pA(x,t) p A (x + f,t) - c A (t) 2 , (2) 

x=l 

where Pa(x) denotes, whether there is an A-atom at the 
surface at site x [pa{x) = 1] or not \pa{x) = 0], and 
CA{t) denotes the surface density of A-atoms. d is the 
dimension of the surface. In our simulations, d was 1 or 2. 
A data collapse of these correlation functions for different 
values of D/F is obtained, if the space coordinates are 
rescaled by £d (see Fig. 3 for d = 1), which shows that 
this is the characteristic distance between the A-clusters. 

The mechanism of the periodic self organization can 
most clearly be seen in the limit rg = E/D — > oo , where 
every S-adatom exchanges with the first exchange part- 
ner A it encounters. When layer t is completed and layer 
t+1 begins to grow, the first _B-atoms deposited are likely 
to exchange with A-atoms from layer t. Hence the nuclei 
of islands in the new layer t+1 will consist predominantly 
of A-atoms. As growth proceeds, the lower terrace (layer 
t) gets depleted from exchange partners A, either because 
they are exchanged with freshly deposited -B-adatoms or 
get overgrown by the islands. Then the core of the islands 
with a high concentration of A-atoms gets surrounded by 
mainly B-atoms (Fig. 2, left). However, as the island 
size increases, it becomes more and more likely that B- 
atoms are deposited on top of the islands, i.e. in layer 
t + 2. These £>-adatoms find many exchange partners 
at the core of the islands, which then become "washed 
out" and start decorating the island edges, because there 
are no Ehrlich-Schwoebel barriers in our model. Note 
that this edge decoration with A-atoms happens with- 
out lateral exchange of B- with A-atoms at the edges of 
the islands, in contrast to the situation studied in [17]. 
As a result, the interior of the islands gets cleared from 
A-atoms, which are collected in the holes of layer t + 1 
which get filled last (Fig. 2, right). Then the process 
starts again: Layer t + 2 nucleates predominantly with 
A-atoms which were exchanged from layer t + 1. 

For finite te the mechanism is similar. However, for 
vicinal surfaces growing in step flow mode, the correla- 
tions among the A-atoms are different. Here, the terraces 
get cleared of A-atoms, which attach to the step edges. 
This decoration of advancing edges leads to a correlation 
pattern g(f,t) with a spatial periodicity identical to the 
width of the terraces. 



IV. THE SCALE FREE LIMIT 

The limit of perfect layer- by-layer growth, D/F — > oo, 
is particularly instructive. In this case, there is only one 




FIG. 2. Top view on the surface structure at t = 3.3 
ML (left) and t = 4.0 ML (right). A-atoms are black, 
B atoms are height-encoded, where brighter means higher. 
D/F = 10 7 , E/D = 10 3 , L 2 = 100 x 100. 
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FIG. 3. Rescaled lateral A-A correlation functions g(r, t) at 
t = 1 for different D/F (d = 1). The average distance between 
A-clusters scales with D/F like the distance of the nucleation sites, 
Its- 



island on the surface, nucleating at a random site. After- 
wards at most one adatom can be found on the surface 
at a time. If we assume maximal exchange in addition, 
i.e. an exchange rate E much larger than both D and 
F, te — E/D — > oo, the model becomes parameter-free. 
A-atoms can be buried even in this limit, when they are 
overgrown by an island edge. However, the last A- atom 
will never get buried in this case. The nucleus of a new 
layer will always contain at least one of the remaining 
surface atoms of type Ain the limit we focus on, because 
the first 5-atom deposited after completion of a layer 
exchanges with an A-atom before the next 5-atom gets 
deposited: ca{P — ► oo) = L~ d . As this case is scale free, 
we expect that the concentration of A-atoms falls off like 
a power law into the growing 5-film. 

In order to get a first idea about the distribution of 
A-atoms in the growing 5-film, it suffices to consider the 
concentration CA(n) of A-atoms at the surface after the 
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deposition of n monolayers. During the growth of the 
next layer, a certain fraction (1 — q) e [0, 1) of these 
A-atoms is transported to the next layer via vertical ex- 
change: 



c A (n + 1) = (1 - q)c A {n). 



(3) 



If q was constant, this would imply an exponential de- 
cay CA(n) oc (1 — q) n . In the present model, however, 
the probability q that an A- atom gets overgrown de- 
creases with decreasing ca(ti), resulting in a decay which 
is slower than exponential. 

Qualitatively this can be understood in the following 
way: In the limit D/F — ► oo there is only one island on 
the surface. Moreover, any A-atom on the lower terrace 
gets transferred to the new layer as soon as it is reached 
by an adatom of type B, because r~E — > oo is assumed, as 
well. Only A-atoms sufficiently close to the island have a 
chance to be overgrown by the island edge before being 
visited by a £?-atom. 

The island edge advances a characteristic distance 
^covcr, while the lower terrace gets depleted from A- 
atoms by exchange with adatoms of type B. For a one- 
dimensional surface of length L it is clear that € C ovor is 
proportional to the number of A-atoms at the surface, 
ca(ti)L. Hence, the number of A-atoms with a chance to 
be overgrown is of the order of 



t ca(ti) oc CA(n) 2 L. 



(4) 



This must be compared with qcA{n)L, which shows that 

qccc A (n). (5) 

Inserting Eq.(5) into Eq.(3) leads to the difference 
equation 



c A (n + 1) - c A (n) oc -c A (nY 
implying the asymptotic power law 
c A {n) ~ 1/n 



(6) 



(7) 



for the concentration of A-atoms at the surface. 

The concentration profile of A-atoms inside the grown 
film is given by 



p(n) = c A (n) - c A (n + 1) ~ 1/n 2 



(8) 



It is remarkable that for the mechanism discussed in this 
paper the width of the interdiffusion zone diverges loga- 
rithmically with the thickness T of the film: 



T-l 

£ 

71=1 



p(n)n + ca(T)T ~ InT. 



(9) 



Nevertheless, the interface can be localized precisely, be- 
cause B-atoms do not occur below layer n = due to the 
absence of bulk diffusion in this model. Below we show 
that these power laws are confirmed by simulations. 



The argument leading to Eq.(5) ignores that the A- 
atoms at the surface are clustered, as shown in Fig. 2, 
and was made plausible only for a one-dimensional sur- 
face. However, it can be refined such that it takes these 
spacial correlations into account and applies also for two- 
dimensional surfaces. For D/F — > oo, we can imagine 
that there is only one cluster of size CA{n)L d on the sur- 
face, when the new layer nucleates. The important point 
is that the nucleation happens anywhere on the surface 
with equal probability l/L d in this case. However, only 
if the nucleation site is within an area of about the size 
2 d CA (n)L d centered at the middle of the cluster, there 
is a chance that some A-atoms get overgrown. In other 
words, only a fraction of nucleation sites oc ca(ti) leads to 
overgrowth. The average number of A-atoms overgrown 
in such a case is proportional to the cluster size. Hence 
on average a fraction q of A-atoms is overgrown which is 
proportional to the fraction of nucleation sites leading to 
overgrowth, i.e. this refined argument gives q oc ca{ti) as 
in Eq.(5). 



V. THE WIDTH OF THE INTERDIFFUSION 
ZONE FOR D/F -> oo 

In the previous section we predicted that the width 
of the interdiffusion zone diverges in the scale free limit, 
where D/F — ► oo and — > oo. In this section we pre- 
dict, that for finite Te the power law Eq.(7) is only valid, 
if the system is infinitely large. For finite system size 
the power law is exponentially cut off at a characteristic 
width of the interdiffusion zone. 

After the completion of several monolayers on a sub- 
strate of linear size L the number of substrate atoms at 
the surface is CAL d . The A-atoms are concentrated in 
a cluster, which we assume to be compact, hence of di- 
ameter oc cJf'L. This assumption is justified even for 
d = 2, where the islands initially are fractal, because the 
A-cluster occupies the sites which were filled last in the 
uppermost monolayer. These sites do not form a fractal. 

Now we imagine the surface to be coarse grained on 
the scale of the cluster diameter so that exactly one cell 
contains the A-cluster. The typical residence time of an 
adatom in such a cell is 



At = 



(c A ,d L? 
D 



(10) 



A -B-adatom which enters the cell containing the A- 
cluster will almost certainly be replaced by an exchange 
partner A within the residence time, if 



EAt = r E (c A /d L) 2 ^>l 



(11) 



(exchange dominates) . For EAt <C 1 the adatom changes 
from type B into type A only with probability EAt (over- 
growth dominates). 

It is plausible to assume that the power law belongs 
to the exchange dominated slow decay of ca while the 
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exponential cut off indicates the much faster decay when 
overgrowth dominates. Thus the width H of the inter- 
diffusion zone should be reached, when ca becomes so 
small that exchange is no longer guaranteed, i.e. when 
EAt drops below 1. Inserting ca ~ 1/t = 1/H into 
Eq.(ll) one obtains 

H^(^¥EL) d . (12) 

However, as ca cannot become smaller than L~ d , this es- 
timate is only valid for Te < 1, while H w L d for Te > 1- 
This is our prediction for the width of the interdiffusion 
zone in the limit D/F — > oo. Note that in this limit the 
cutoff of the power law is a finite size effect: For L —* oo 
the power law extends to infinity. 

Based on the results of this and the previous para- 
graph we can conjecture the following scaling form for 
the surface concentration of ^4's : 

c A (t,L;r E )-c A (t^^) = ^f(j^j (13) 

where according to Eq.(7) 

/(t)~1/t for t<1. (14) 



and varied the values of te and system size L for one- and 
two-dimensional surfaces. For the case d = 1 we used 
the algorithm described in the appendix, while kinetic 
Monte-Carlo simulations were done for d = 2. 

Fig. 4 (for d = 1) and Figures 5, 6 (for d = 2) show the 
concentration ca of ^4-atoms at the surface as a funcion of 
deposition time t (in monolayers of 5-atoms) . All curves 
are averages over 200 - 400 independent runs. Both for 
d = 1 and d = 2 the exponent of the power law decay 
was found to be consistent with the value -1 derived in 
section IV. 

As shown in Fig. 4, the predicted relations Eq.(12) 
and Eq.(13) lead to the expected data collapse for the 
one dimensional surface. The results in two dimensions 
are not as clear. In this case, we obtain the best data 
collapse with 

H oc L 1 ' 93 • 4' 2 , (15) 

as shown in the Figures 5 and 6, whereas our predicted 
exponents (2 and 1, respectively, see Eq.(12)) were about 
4% and 20% different. In fact the ^4-clusters are not as 
compact as assumed in the simple argument of Sec.V (see 
Fig.2). 
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t / (L • (E/D) 5 ) 

FIG. 4. Scaling of c A (t) with L and E/D in d = 1 for 
D/F — > oo. Full symbols mark four curves for E/D = 0.1 and L 
between 256 and 8192. Open symbols mark five curves for L = 4096 
and E/D between 10 — 4 and 100. A data collapse of all 9 curves 
is reached by scaling in accordance with Eq.(13) and Eq.(12). The 
data for E/D = 100 were rescaled differently: Here, Cj\L — 1 is 
plotted vs. t/L, as if E/D was 1 instead of 100. This shows, that 
H becomes independent of E/D for E/D > 1 as explained after 
Eq.(12). The dashed line has slope -1 in agreement with Eq.(14). 



VI. NUMERICAL RESULTS FOR D/F — > oo 

In order to check the predictions of Sees. IV and V, we 
simulated the model described in Sec. II for D/F — > oo 
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FIG. 5. 

D/F — > oo. 
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Scaling of CA(t) with te 
L 2 = 200 X 200. 



E/D in d = 2 for 



VII. SCALING FOR FINITE D/F 

For finite D/F there are many A-clusters on the sur- 
face. Their typical distance is given by the diffusion 
length £u as shown in Sec. III. The average size of the A- 
clustcrs is ca^d, provided this is much larger than 1. If 
the concentration ca becomes too small, less and less A- 
clusters will be found on the surface, and their typical dis- 
tance will grow. Finally all ^4-atoms will be overgrown, in 
contrast to the situation of perfect layer-by-layer growth, 
where the last A-atom could never be overgrown. Apart 
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FIG. 6. Scaling of c A (t) with L in d = 2 for D/F 
E/D = 0.1. 
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FIG. 7. 

E/D = 10 3 ,L = 5 
exponent j3 = 0.78. 



Scaling of ca(<) with D/F in d = 1. 

10 4 . The dashed line indicates the 
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from this, one might expect that the results of the pre- 
vious three sections would essentially remain true, if one 
replaces L by in- Qualitatively, the surface concentra- 
tion ca indeed decays first approximately as a power law 
of the deposition time, which is cut off at a characteristic 
width H of the interdiffusion zone. 

Quantitatively, however, the situation turns out to be 
more complex than this: All exponents are different, as 
the simulation results (Figs. 7 - 10) show. The power 
law decay of ca oc £ -/3 extends over at most two decades 
for the largest values of D/F we simulated, so that the 
determination of the exponent (3 from the slopes in the 
log-log-plots Fig. 7 and 8 is not very accurate. We esti- 
mate 

f 0.78 ±0.08 for d=l, , . 

\ 0.53 ±0.05 for d = 2, ! ' 

which arc indicated by the dashed lines in the two Fig- 
ures. Both exponents are significantly smaller than (3=1 
obtained for infinite D/F, i.e. ca decays more slowly for 
finite than for infinite D/F. 

This result is surprising on first sight, because there are 
more island edges on the surface for finite D/F, hence 
more possible places, where ^4-atoms may be overgrown. 
That ca decays more slowly nevertheless, may be ex- 
plained by the fact, that the nucleation of islands does 
not happen anywhere with equal probability as for infi- 
nite D/F , but preferentially far away from the holes in 
the previous layer, where the ^4-atoms are concentrated. 
Therefore overgrowth is less likely, and ca decays more 
slowly than for infinite D/F. 

This raises the question, how big the parameter D/F 
must be in order to see the exponent (3=1 instead of 
the smaller one. The answer is, that the system size 
L must be small compared to £n in order to obtain the 
crossover to the faster decay of ca- This was confirmed by 
simulation results in [18]. In the four Figures belonging 



to this Section we carefully checked that the system sizes 
were big enough to exclude finite size effects. 

In analogy to Eq.(ll) we expect that the power law 
decay of ca oc t _/3 stops, when 

Mc A ^f /d « 1, (17) 

or 

CAlfy « 1, (18) 

whichever happens first. Replacing ca by , this im- 
plies that the width H of the interdiffusion zone should 
be given by 

H^£ d / j for r E »l (19) 

and 

H^(^tu) d/I3 for r E «l. (20) 
In analogy to Eq.(13) we postulate then that 

- = h* (1) ^ 

with 

9 {t) - l/r f3 for t < 1, (22) 

because ca is independent of H for small t. 

We first checked these conjectures for r E > 1, where 
the surface concentration of A 1 a becomes independent of 
r E , as expected. If we insert the D/i^-dependence Eq.(l) 
of the diffusion length in Eq.(19), Eq.(21) can be written 
in the form 

-(§r=4(fn- 
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FIG. 8. 

E/D = 10 4 ,L 2 
poncnt (3 = 0.53. 
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Scaling of c^(t) with D/F in d = 2. 
= 500 X 500. The dashed line indicates the ex- 
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With 7=1/4 [19] for d = 1 and 7 = l/(4 + d f ) w 0.17 for 
d = 2 (df is the fractal dimension of the islands) [20] , and 
with the /3-valucs determined above, the theory predicts 



7c? = 



id/ /3 = 



0.25 for d=l 

0.35 ±0.01 for d = 2 



0.32 ± 0.03 for d=l 
0.66 ± 0.07 for d = 2 



(24) 



(25) 



The data collapses in Figs. 7, 8 are in reasonable agree- 
ment with this prediction. 

However, the rE-dependence Eq.(20) for te < 1 is not 
in agreement with the simulation results. For fixed D/F 
the theory Eq.(21) predicts 



d/2 

c A r^ 



-d/20s 



(26) 



Inserting the value of (3 determined above, the scaling 
exponent should be 



d/20-. 



0.64 ±0.06 
1.9 ±0.2 



for d = 1 
for d = 2 



(27) 



For D/F — 10 7 we could only check this for about one 
decade of rs-values: For d = 1 we found that already for 
te = 0.3 the crossover into the regime, where H becomes 
independent of r&, affects the data. For te < 10~ 3 the 
exchange was so weak that the surface concentration of 
j4's decayed very fast from the beginning, so that a con- 
vincing data collapse was not possible. Similar problems 
occurred for d=2. The best result of our attempts to get 
a data collapse in the available rE-interval are shown in 
Fig. 9 for d = 1 and Fig. 10 for d — 2. The effective 
exponents turn out to be very different from the ones 
predicted in Eq.(26). 
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VIII. RATE EQUATION APPROACH 

In this Section we extend the established rate equation 
approach for submonolayer homoepitaxial growth as de- 
scribed in [21,22], in order to apply it to our model for 
surface interdiffusion. Our approach describes the time 
evolution of four submonolayer quantities: The density 
of mobile adatoms p, the total island density /, the den- 
sity of mobile adatoms of type B, pb, and the density of 
potential exchange partners of type A in the lower layer, 
Pa- With these quantities, the rate equations are as fol- 
lows: 

p = F-Dp(I + 2p), (28) 
I = Dp 2 , (29) 
p B = F- Dp B (I + p + p B )- Ea 2 p A p B , (30) 
p A = -Ea 2 p A p B - p A a 2 Dp{I + 2p). (31) 
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While the first two equations are identical to the well 
known point-island model rate equations for homoepitax- 
ial growth, the last two are specific to our hctcrocpitaxial 
model. The third one expresses the change in density of 
the -B-type adatoms. Its positive contribution describes 
the depostition of new adatoms. The first negative term 
represents the loss of 5-adatoms, when they get incor- 
porated into islands or bind to another adatom to nu- 
cleate a new island. The extra term Dp 2 B accounts for 
the fact that nucleation events involving two B-adatoms 
count twice as much as those between a B- and an A- 
adatom, because they remove two 5-adatoms simultane- 
ously. The last term of Eq.(30) describes the exchange of 
mobile B's with A's. Equation (31) expresses the annihi- 
lation of possible A-type exchange partners in the lower 
layer. Since this value is monotoneously decreasing, there 
is no positive contribution. The negative terms describe 
the exchange of A's with mobile B's, and the overgrowth 
of A's due to propagating island edges and nucleation 
events. 

Rescaling the variables (see [21]) according to i 

p = pl\, I = lt\, pA = PA?-l, PB = Pb^o, where £ 
(D/F) 1 / 4 , leads to the dimensionless equations 



tFl%, 



p=l-p(I + 2p), 
I 

PB 



P 
1 



Pb(I + f> + Pb) -r E pAPB, 

2 

PA = -r E pAPB - (a/l ) f>Ap(I + 2/3). 



(32) 

(33) 
(34) 
(35) 



If we consider systems in perfect layer-by-layer growth 
mode, this approach not only holds for the submonolayer 
regime starting from the substrate, but also starting af- 
ter integer numbers of deposited monolayers. The initial 
conditions of these equations for a flat surface after n 
deposited monolayers are 

p(0) = 1(0) = /3 B (0) = 0, p A (0) = c A (n). (36) 

Disregarding the point island model nature of this ap- 
proach, which only holds for early stages of the submono- 
layer regime, we can establish the surface concentration 
of ^4's after the deposition of one additional monolayer, 
CA(n+l), as the integral over the density of all exchanged 
atoms: 



c A (n + 1) 



-I 



(to/a) 2 



di r E p A Pb 



(37) 



The upper integration boundary is the dimensionless 
time for depositing one monolayer. This approximation 
can be justified by taking into account that the transport 
of ^4-atoms from the nth to the (n + l)th layer mainly 
takes place at early times, that is, the nucleation regime 
and early stages of the intermediate coverage regime, as 
explained in section III. The chosen approach describes 
these regimes with sufficient accuracy. 

The solution of the first two equations can be taken 
directly from the literature [22]: For early times, t< 1, 



p is linear in i, and / increases with P. At late times, 
i » 1, one gets p cx i 1 / 3 and / cx £ -1 / 3 . 

With these results, the last two equations can be solved 
analytically in a similar way for the early-time regime, 
( « 1. For equation (34), the second term on the right 
hand side can be neglected in this limit. If we also neglect 
the time dependence of pA{t) ~ Pa{0) = CA(n), we get 



Pb ~ 1 - r E ca Pb 



(38) 



This equation relaxes into a steady state with ps,oo °^ 
V(rE ca) after a characteristic time t* ~ l/(r E CA(n)). 
For even smaller times, t C f, we can also neglect the 
other right hand side term, and we get p E oc t. 

Plugging these results into Eq.(35), and realizing that 
the second term on the right hand side can be neglected 
compared to the first one, we get 



PA oc e 

and 

PA oc e 



-i/c A 



-r E t A 



(1 - r E P) 



for 



i* < t < 1 (39) 



for 



(40) 



To relate these findings to our results in the other 
sections, we employed an iteration scheme for the rate 
equation system to obtain the surface concentration of 
A's for every deposited integer monolayer. Starting 
from the substrate (pa(0) = ca(0) = 1) we can obtain 
Pa{1) = c A (l) from Eq.(37) by solving Eqs.(32) - (35) 
numerically. Plugging ca(1) back into our rate equations 
as the initial surface concentration, that is pa(0) — ca(1), 
we get pa(1) = c^(2) by using Eq.(37) again. Repeating 
this iteration scheme leads to CA{t). 
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FIG. 11. Scaled time dependence of the surface concentration 
of substrate atoms, obtained from the iteration of the rate equa- 
tions. The slope of the power law region is —1.03 ± 0.05. E/D = 
1. 



The rescaled results of this approach for different D/F 
are shown in Fig. 11. One can clearly observe power law 
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decay for high D/F values at intermediate times, and 
a similar scaling behaviour as obtained from the sim- 
ulations. The exponent of the power law behaviour is 
approximately —1, which is identical to the result from 
the simulations for perfect layer- by-layer growth mode, 
D/F — > oo. This fact supports our argumentation con- 
cerning the different exponents for D/F — > oo and finite 
D/F in section VII: Since the rate equations cannot de- 
scribe the clustering of the yl-atoms, they also do not 
reflect the preference of the nuclcation sites to be far 
from the yl-clusters. The observed scaling exponent of 
Zq = (D/F) 1 ^ 2 doesn't match any of the exponents re- 
sulting from the simulations. This is in general agree- 
ment with the analytical calculation of the exponents in 
section V: There we derived the scaling exponents from 
characteristic properties of the yl-clusters, which are to- 
tally neglected in the presented rate equation approach. 

IX. CONCLUSION 

In the present work we have investigated heteroepitax- 
ial growth of 5-particles on an ^4-substrate. Introducing 
an exchange mechanism for B-adatoms, when they en- 
counter an ^4-atom in the uppermost layer, we observed 
that in the limit of layer-by-layer growth the top layer 
concentration of yl-atoms decays algebraically. There- 
fore, the resulting interdiffusion zone has a broad profile 
with a diverging width. Varying the rate E at which 
^4-atoms and 5-atoms are exchanged did not change the 
exponent of this power law. A different situation has 
been observed, as we varied the diffusion constant D. 
For finite values of D a crossover from power law to ex- 
ponential decay has been found. The crossover time H is 
given by (D/F) - 3 for d = 1 and by (D/F) 6 for d = 2. 

It would be very interesting to have some experimental 
results for the concentration of substrate atoms in the 
growing layer in the case of heteroepitaxy. The results we 
found suggest that in order to provide a sharp interface 
between substrate and deposited material one should use 
small D/F. 
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APPENDIX 

Here we describe, how we implemented the model in- 
troduced in Sec. II for one dimensional surfaces in the 
limit D/F — > oo. We first describe the idea for the scale 
free limit, where also E/D — > oo. 

For D/F — > oo one has perfect layer-by-layer growth. 
The nucleation of a new layer happens at an arbitrary 
position. Afterwards there is at most one adatom on the 
surface. The idea is to calculate the probabilities exactly, 
with which the adatom reaches the nearest sinks to its 
left and to its right. For an ^4-adatom these are the island 
edges, while for a B-adatom it might also be an yl-atom, 
with which it could exchange. Let dh (dn) denote the 
distance to the nearest sink to the left (right). 

As shown in [23], the probability pl to reach the left 
position prior to the right one with unbiased diffusion is 
given by 



Correspondingly, p R = 1 — p L - Therefore, it is not nec- 
essary to simulate the whole random walk of an adatom, 
but it suffices to select the final position according to 
(41). 

Thus in the scale free limit the model (after nucleation 
of a new layer) may be simulated as follows: 

(1) Deposition of a Bat a randomly chosen site i. 

(2) Determination of the distances dh and cZr followed by 
a decision for a side according to the probabilities given 
in (41). 

(3) If the final position of the B-adatom is an ^4-site, the 
atoms exchange (as E/D — > oo). In this case the A- 
adatom goes to the left or right island edge according to 
(41). If the final postition of the B-adatom is an island 
edge, it is bound there irreversibly possibly overgrowing 
an v4-atom. Then one returns to step (1) and deposits 
the next B-atom. 

This algorithm can be generalized for finite te = E/D: 
Not always, when a B-adatom encounters an exchange 
partner A, they exchange immediately. This happens 
only with probability pe = E / (E + 2D), where the de- 
nominator is the sum of the rates for the three possible 
actions of the adatom - exchange with the A-atom under- 
neath, a hop to the right neighbour and a hop to the left 
neighbour. With probability pe the -B-adatom is replaced 
by an ^4-adatom, which attaches to the island edges to its 
left with probability Eq.(41), and otherwise to the island 
edge to its right; 1 — pe is the probability that the B- 
adatom continues to diffuse until it encounters the next 
^4-atom or attaches to the island edge. 

In order to avoid simulating the random walk ex- 
plicitely, one has to calculate the probabilities analyti- 
cally, that the B-adatom exchanges with any particular 
of the ^4-atoms or attaches to the island edges. Techni- 
cally speaking, the B-atom is a random walker on a one- 
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dimensional lattice with fixed partial absorbers (the A- 
atoms) and two full absorbers (the island edges) (Rosen- 
stock trapping model with partial absorbers). In order to 
calculate the absorption probabilities at the different ab- 
sorbers, which depend on the deposition site, we consider 
an incoming flux (normalized to 1) of independent ran- 
dom walkers at the deposition site xs (source) and deter- 
mine the outgoing fluxes at the absorption sites (sinks). 
The absorption probability is then the steady state frac- 
tion of the incoming flux that leaves the system at the 
respective absorption site. 

The density of random walkers at a site x evolves ac- 
cording to 

p(x, t) = D[p(x - 1, t) - 2p(x, t) + p(x + 1, *)] (42) 
-Ep(x,t)p A (x) + d x<xs , 

where the density of partial absorbers, Pa(x), is 1 at all 
the sites x\, where an yl-atom sits, and otherwise: 

n 

Pa(x) = ^^, IA „. (43) 

The terms on the right of Eq.(42) which are proportional 
to D are the gain and loss terms due to hopping from a 
neighbor site to x, respectively away from x. The term 
proportional to the exchange rate E describes the loss 
of walkers at the partial absorption sites. The last term 
is the gain term due to the normalized influx of walkers 
at site xs ■ The perfect sinks corresponding to the island 
edges are represented by the boundary conditions p(l) = 
p(L) = 0, where L is the size of the terrace, on which the 
source is located. 

The probability of absorption at site xa is then ob- 
tained from the steady state solution of Eq.(42) by 

p(x A ) = Ep(x A ), (44) 

and the ones at the island edges by 

p(l) = £>p(2), p(L) = Dp(L-l). (45) 

Introducing the diffusion current between x and x + 1 
(i.e. the current to the right of x and to the left of x+ 1), 

jn(x) = 3l(x + 1) = -D(p(x + 1) - p(x)), (46) 

equation (42) can be rewritten in the steady state as 

n 

Jr(x) - jh(x) = -Ep{x) ^2 ^x.xa, + £x,x s - (47) 

v=\ 

This shows that p(x) is a piecewise linear function with 
slope discontinuities at the source and the sinks. Hence 
Eq.(42) reduces to a set of 2n + 2 coupled linear equa- 
tions for the 2n + 2 unknowns .)r(xaj/), p{xau) and the 
boundary values Jr(1) and jh(L). 

The solution determines the probabilities Eqs.(44), 
(45) with which a freshly deposited 5-atom is exchanged 



at the different ^4-sites or absorbed by the island edges. 
By choosing a random number we decide which site to 
pick. If it is an island edge, the 5-atom is moved there, 
and the next 5-atom is deposited at a random position. 
Otherwise we move the 5-atom to the chosen site, ex- 
change it with the yl-atom there, let another random 
number determine, whether to attach the ^4-atom to the 
left or right island boundary, and deposit the next B- 
atom at a random position. 

The complexity of this algorithm is linear in the num- 
ber of ^4-atoms left on the surface, while a brute force 
simulation of the diffusion would cost much more com- 
puting time proportional to L 2 . 
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